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On the X-Ray and Mass Distribution in the Merging Galaxy Cluster IE 
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Abstract 

Wc investigate the X-ray and mass distribution in the merging galaxy cluster IE 0657-56. We study 
head-on collisions of two virialized clusters with an NFW density profile in the ACDM universe using an 
TV-body + hydrodynamical code. A clear off-set of an X-ray peak from a mass peak, which is like what 
is reported in Clowe et al. (2004), is first reproduced in the A"-body + hydrodynamical simulations. We 
estimate the ram pressure-stripping conditions of the substructure in mergers of two NFW dark halos 
using a simple analytical model. We find that the ram pressure dominates the gravity of the substructure 
when the smaller cluster's mass is less than approximately one tenth of the larger cluster's mass. The 
characteristic X-ray and mass structures found in IE 0657-56 suggest that neither the ram pressure nor 
the gravitational bound force overwhelms the other and that the mass ratio between the progenitors is 
near the critical value mentioned above. 
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1. Introduction 

According to the standard scenario of cosmological 
structure formation, clusters of galaxies form through 
mergers of smaller subclustcrs. In fact, some of them are 
forming now. Merging clusters are the sites of structure 
formation in the universe that can be investigated in detail 
via different types of observations. Mergers play impor- 
tant roles in cluster evolution itself. Cluster mergers cause 
bulk flow motion, turbulence, shocks, and/or contact dis- 
continuities in the intraclustcr medium (ICM) , which give 
us clues to investigate various physical processes in the 
ICM (see Sarazin 2002 for a review). Turbulence and 
shocks most likely play crucial roles in particle acceler- 
ation in the ICM (e.g. Takizawa, Naito 2000; Ohno, 
Takizawa, Shibata 2002; Takizawa 2002; Fujita, Takizawa, 
Sarazin 2003; Brunetti et al. 2003 ). Strong bulk flow 
motion and violent pressure changes in the ICM during 
mergers may affect star formation activities of the member 
galaxies (Fujita et al. 1999; Bekki, Couch 2003). N-body 
+ hydrodynamical numerical simulations have been car- 
ried out to study cluster merger physics (e.g. Schindler, 
Miiller 1993; Ishizaka, Mineshige 1996; Roettiger, Loken, 
Burns 1997; Takizawa 1999, 2000; Ricker, Sarazin 2001; 
Ritchie, Thomas 2002). Comparison of such numerical 
simulations with different kinds of observations give us 
deep insight into the cluster physics. 

IE 0657-56 is one of the most well-known examples of 
a merging cluster. It is the hottest known cluster and has 
a very powerful radio halos (Liang et al. 2000). Its sim- 
ple geometry makes this cluster one of the most suitable 



case to investigate cluster merger physics. There are two 
peaks in both the X-ray surface brightness distribution 
(Markevitch et al. 2002) and galaxy distribution (Barrena 
et al. 2002), but their positions do not agree with each 
other. Observations of the line-of-sight velocities of the 
member galaxies suggest that its collision axis is almost 
perpendicular to the line-of-sight (Barrena et al. 2002). 
From X-ray observations, a bow shock and a cold front 
are found in front of the smaller subcluster, which suggest 
that its Mach number is 2 — 3 (Markevitch et al. 2002). 

Recently, Clowe et al. (2004) investigated the mass dis- 
tribution in IE 0657-56 through weak gravitational lens- 
ing. They show clear offsets of the mass density peaks 
from the X-ray peaks, and that the mass distribution is 
quite similar to the galaxy one. The smaller substruc- 
ture in mass is ahead of the X-ray one. They claim that 
this structure occurs because ICM experience ram pres- 
sure but dark matter (DM) and galaxies do not. Although 
the above-mentioned naive ram pressure-stripping sce- 
nario seems to be correct, such characteristic off-sets of 
X-ray peaks to mass peaks have never been reported in 
the past numerical simulations. In this paper, wc show the 
first results that successfully reproduce such characteristic 
structures in the Af-body + hydrodynamical simulations, 
and discuss their implications using a simple analytical 
model. 

The rest of this paper is organized as follows. In §2 wc 
describe the adopted numerical method and initial condi- 
tions for our simulations. In §3 we present the simulation 
results. In §4 we show simple analytical estimation for 
ram pressure-stripping conditions in mergers of two clus- 
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ters with an NFW density profile. In §5 we summarize 
the results. 

2. The Simulations 

2.1. The Numerical Method 

In the present study, we consider clusters of galaxies 
consisting of two components: collisionless particles cor- 
responding to the galaxies and DM, and gas correspond- 
ing to the ICM. When calculating gravity, both com- 
ponents are considered, although the former dominates 
over the latter. Radiative cooling and heat conduction 
are not included. We use the Roe total variation dimin- 
ishing (TVD) scheme to solve the hydrodynamical equa- 
tions for the ICM (see Hirsch 1990). The hydrodynam- 
ical part of the code used here is identical with what is 
used in Takizawa (2005). Gravitational forces are calcu- 
lated by the Particle-Mesh (PM) method with the stan- 
dard FFT technique for the isolated boundary conditions 
(sec Hockncy, Eastwood 1988). The size of the simulation 
box and the number of the grid points arc 11.8 Mpc x 
(5.92 Mpc) 2 and 256 x (126) 2 , respectively. The total 
number of the iV-body particles used in the simulations is 
256 x (128) 2 , which is approximately 4.2 x 10 6 . 

2.2. Models and Initial Conditions 

We consider mergers of two virialized subclusters with 
an NFW density profile (Navarro, Frenk, White 1997) 
in the ACDM universe (fi = 0.25, A = 0.75) for DM. 
DM masses of the larger and smaller subclusters are 
1.00 x 1O 15 M and 6.25 x 1O 13 M , respectively. Thus, 
the mass ratio is 16:1. The initial density profiles of the 
ICM are assumed to be those of a beta-model, where the 
core radius is a half of the scale radius of the DM distri- 
bution, and (3 = 0.6. The gas mass fraction is set to be 0.1 
inside the virial radius of each subcluster. The parame- 
ters such as a virial radius r v , concentration parameter c 
for each subcluster are summarized in table 1. We calcu- 
late these parameters following the method in Appendix 
of Navarro, Frenk, White (1997). The radial profiles of the 
ICM pressure is determined so that the ICM is in hydro- 
static equilibrium within the cluster potential of the DM 
and ICM itself. The resultant temperature profiles are 
similar to those of initial cluster model in Ricker, Sarazin 
(2001). The velocity distribution of the DM particles is 
assumed to be an isotropic Maxwellian. The radial pro- 
files of the DM velocity dispersion are calculated from the 
Jeans equation with spherical symmetry so that the DM 
particles are in virial equilibrium in the cluster potential. 
The coordinate system is taken in such a way that the 
center of mass is at rest at the origin and that the x— axis 
is along with the collision axis. The centers of the larger 
and smaller subclusters are initially located in the sides of 
x > and x < 0, respectively. The initial distance between 
the subcluster's centers is 4.93 Mpc. The initial relative 
velocity is estimated as in §2 of Takizawa (1999). The re- 
sultant value is 8.98 x 10 2 km s _1 , which is approximately 
two thirds of the infall velocity assuming that they were 
at rest at infinite distance. 



Table 1. Parameters for each subcluster 





M DM (M q ) 


r v (Mpc) 


c 


cluster 1 


1.00 x 10 ib 


1.97 


5.66 


cluster 2 


6.25 x 10 13 


0.784 


7.56 



3. The Simulation Results 

Figure 1(a) shows the X-ray surface brightness distri- 
bution (colors) and projected total surface mass density 
(contours) at a time of 0.67 Gyr after the passage of the 
subcluster through the core of the larger one. A clear 
off-set of the mass density peak from the X-ray peak is 
seen for the smaller subcluster remnant. The mass peak 
and X-ray peak are located at x ~ 1.5 Mpc and x ~ 1.0 
Mpc, respectively. This is clearly because the ICM in the 
smaller subcluster is lagged by the ram pressure. Figure 
1(b) shows the X-ray surface brightness distribution (con- 
tours) overlaid with the emissivity-wcightcd temperature 
distribution (colors) at the same epoch. A weak jump in 
the X-ray surface brightness distribution at x ~ 1.5 (near 
the smaller mass peak) is a bow shock, and the emissivity- 
wcightcd temperature is higher in the brighter side, and 
vice versa. A more prominent jump in the X-ray bright- 
ness just in front of the smaller X-ray peak is a contact 
discontinuity, and the emissivity-weighted temperature is 
lower in the brighter side, and vice versa. Therefore, this 
will be recognized as a cold front in actual X-ray obser- 
vations. Figure 2(a) and (b) show the ICM density and 
pressure profiles along the collision axis (y = z = 0) in front 
of the smaller X-ray peak, respectively. The bow shock is 
located at x ~ 1.55 Mpc, where jumps are clearly seen in 
both the density and pressure profiles. The contact dis- 
continuity is at x ~ 1.1 Mpc, where a jump is seen only in 
the density profile and the pressure profile does not have 
any discontinuity there. As for the overall ICM and DM 
structures of IE 0657-56 around the west smaller X-ray 
and mass peak, our results agree qualitatively with the 
observations. 

4. Discussion on the Ram Pressure-Stripping 
Conditions 

Let us discuss the ram pressure-stripping conditions in 
merger of two clusters with an NFW DM density pro- 
file. We consider the merger of two clusters with masses 
Mi and (M\ > M2), respectively. We concentrate 
on the physical status of the ICM in the core region of 
the smaller subcluster. If the gravity on the subcluster's 
ICM is weaker than the ram pressure force in unit volume, 
the ICM will be stripped from the substructure potential. 
This means, 

^<A(nrl Pl v*)(inrl)-\ (1) 
"2 

where G is the gravitational constant, and p\ and P2 are 
the central gas density of the subcluster 1 and 2, respec- 
tively. r 2 and mi are the scale radius of the DM profile 
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Fig. 1. (a)Projected total surface mass density (contours) 
overlaid with X-ray surface brightness distribution (colors) 
at a time of 0.67 Gyr after the core passage. A clear off-set of 
the mass density peak from the X-ray peak is seen. (b)X-ray 
surface brightness distribution (contours) overlaid with emis- 
sivity-weighted temperature distribution (colors) at the same 
epoch. 
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Fig. 2. Electron number density (a) and pressure (b) profiles 
along the collision axis (y = z = 0) in front of the substructure 
at the same epoch as in figure 1. A bow shock is located 
at x ~ 1.55 Mpc, where jumps are clearly seen in both the 
density and pressure profiles. A contact discontinuity is at 
x ~ l.lMpc, where a jump is seen only in the density profiles 
and the pressure profile has no discontinuity there. 

and the DM mass inside r 2 for the cluster 2, respectively. 
Therefore, the relation between m 2 and M 2 is 



g{M 2 ) 



m 2 
M2 



ln2-l/2 



(2) 



ln(l + c) - c/(l + c) ' 

where, c is a concentration parameter of an NFW profile, 
and weakly depends on the halo mass. A is a fudge factor 
of an order of unity. It is most likely that A<1 because 
all of the ram pressure force is not effective in stripping 
the gas from the substructure. Some might be used on the 
excitation of small-scale eddies through Kelvin- Hclmholtz 
instability, and some on the adiabatic compression and 
shock heating of the ICM, and so on. The collision velocity 
v has an order of 

^2 G (M 1 + M 2)) (3) 
Ri + R 2 

where R± and R% are the virial radii for the cluster 1 
and 2, respectively. It is convenient to introduce a new 
parameter a = MijM\. Then, the stripping condition of 
inequality (1) becomes 

F(a ]Ml)s a V^l±^L _ < 0(4) 



1 



- Q 



2g(aM 1 )c{aM 1 ) 



Deriving inequality (4), we use the scaling relation that 
R 2 /R 1 = a 1/3 and p 2 /pi = a~ w . In the ACDM universe 
(fio = 0.25, Ao = 0.75), w ~ 1/4 assuming that p\ and 
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Fig. 3. Function F(a : Mi) defined by equation (4) for 
Mi = 1.0 X 1O 15 M in the ACDM universe (Q = 0.25, 
Ao = 0.75). The solid, dashed, and dot-dashed lines repre- 
sent the cases of A = 0.6, 0.4, and 0.2, respectively. F < 
means that the ram-pressure dominates the gravity in the 
substructure. 

P2 behave like a characteristic density in an NFW profile 
(Navarro, Frenk, White 1997). 

Figure 3 shows the function F{a : Mi) for Mi = 1.0 x 
10 15 M Q in the ACDM universe (Q = 0.25, A = 0.75). The 
solid, dashed, and dot-dashed lines represent the cases 
of A = 0.6, 0.4, and 0.2, respectively. In any case, a is 
less than ~ 0.1 when F(a) < 0. This means that the 
ram pressure-stripping is more effective for smaller sub- 
clusters. Although our estimation here is rather crude, 
it is interesting that this criteria of a ~ 0.1 is close to 
the mass ratio of our simulations where the clear off-set 
appears. Obviously, such an off-set does not appear at 
all if the ram pressure-stripping does not work effectively. 
The X-ray peak will correspond with the mass peak be- 
cause the ICM behaves like the DM. On the other hand, 
if the ram pressure overwhelms the gravity, the ICM in 
the substructure will not be able to penetrate the larger 
cluster's center and will be repelled. Furthermore, the 
larger cluster's ICM is so hot that it cannot be bound by 
the substructure's gravitational potential. Therefore, we 
will see a mass peak associated with no X-ray peak. Clear 
off-sets in the simulation and IE 0657-56 suggest that the 
parameter a is close to the above-mentioned critical value 
where neither the ram pressure nor the gravity absolutely 
dominates the other. 

In our simulations, we do not take into account a dense 
cool core in the central region of a cluster. However, 
this probably affects the ram pressure stripping condi- 
tions. The existence of the cool core changes p\ and P2 
in inequality (1). Thus, the scaling relation P2I P\ = ct~ w 
could be modified and/or have some dispersion. In case 
of 1E0657-56, the smaller subcluster certainly has a cool 
core, but the larger one seems to have none. Therefore, 
P2I P\ can be larger than when neither subcluster has a 



cooling core. This makes more difficult to strip the ICM. 

Whereas the smaller X-ray peak in 1E0657-56 is tri- 
angle, that in our simulation results is rather round. The 
absence of a cool core might cause this difference. Another 
physical process that may play a crucial role is the mag- 
netic field in the ICM. The magnetic field along the bound- 
ary layer probably works to maintain the smaller subclus- 
ter's gas as a distinct structure after it is stripped off the 
DM potential through both the magnetic tension and sup- 
pression of heat conduction (Asai, Fukuda, Matsumoto 
2004). Dynamical motion of the substructure itself pos- 
sibly produces this kind of magnetic field configurations 
(Vikhlinin, Markevitch, Murray 2001). Furthermore, tem- 
perature gradients in the boundary layer might produce 
the magnetic field structure through Weibel-type plasma 
instabilities (Okabe, Hattori 2003). Three-dimensional 
high-resolution magnctohydrodynamic simulations will be 
useful in order to investigate detailed evolution. 

Please note that an off-set of an X-ray peak to a mass 
peak is not a structure in dynamical equilibrium but 
a transient one. A characteristic timescale of the ram 
pressure-stripping is estimated to be ~ R2/V, which be- 
comes 0.42 Gyr for the model presented in §2. In the 
simulation, it takes ~ 0.7 Gyr to strip the ICM from the 
substructure potential. Certainly this is roughly equal to 
the above-mentioned value. 

5. Summary 

We investigate the X-ray and mass structures in the 
merging galaxy cluster IE 0657-56. We first reproduce a 
clear off-set of an X-ray peak to a mass peak in A-body + 
hydrodynamical simulations of mergers of two subclusters 
with an NFW and beta-model density profiles for DM and 
ICM, respectively. As for the overall ICM and DM struc- 
tures of IE 0657-56 around the west smaller X-ray and 
mass peak, our simulation results agree qualitatively with 
the observations. We discuss the ram pressure-stripping 
conditions in the mergers of two clusters with an NFW 
DM density profile using a simple analytical model. We 
find that the ram pressure dominates the gravity of the 
substructure when the smaller cluster's mass is less than 
approximately one tenth of the larger cluster's mass. The 
characteristic X-ray and DM structures found in IE 0657- 
56 suggest that the mass ratio between the progenitors 
is close to the above-mentioned critical value and that 
neither the ram-pressure force nor the gravity of the sub- 
structure overwhelms the other. 
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